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Imaging 



5 This invention relates to imaging and in particular to 
methods and systems for image enhancement. 

Imaging involves transfer from the object domain into the 
image domain, but owing to limiting factors such as the 

10 finite size of energy source, detector size, sampling 

frequency, display density, software filter function, and 
possibly partial-volume effects experienced with some 
imagers, an infinitely fine delta function in the object 
domain cannot be faithfully reproduced in the image 

15 domain. Instead, a smeared-out image, or point-spread 
function (PSF) , is observed. Similarly, an infinitely 
sharp edge-response function (ERF) in the object domain 
becomes a smeared-out ERF in the image domain. The 
smearing effect becomes more intense as the adjacent ERFs 

20 of discontinuities or contrast profiles get closer to 
each other. 

It is an object of the present invention to provide a 
method and system by which the above problem can be at 
25 least partly overcome. 

According to one aspect of the present invention there is 
provided a method wherein a de-convolution process is 
applied to the image-domain results of an object-scan to 
30 derive therefrom the respective point- or line-spread 
function effective in the object- to image-domain 
transfer of one or more object-discontinuities, and to 
derive from said function the location in the image 
domain of the respective discontinuity. 

35 

According to another aspect of the invention there is 
provided an imaging system comprising means for 
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performing a de-convolution process on the image-domain 
results of an object-scan to derive therefrom the 
respective point- or line-spread function effective in 
the object- to image-domain transfer of one or more 
object-discontinuities, and means to derive from said 
function the location in the image domain of the 
respective discontinuity. 

The method and system of the invention enable the 
location of the respective discontinuity in the image 
domain, to be established with a high degree of accuracy. 
This is critical to image definition free of any 
substantial smearing, and to this end the location of the 
respective discontinuity may be derived to sub-pixel 
accuracy simply from the mid-point of the full-width 
half -maximum of said function. 

The said function may be correlated with the image-domain 
results of said transfer for enhancement of spatial 
resolution of the imaging of the one or more 
discontinuities. This enhancement may involve transfer 
of sub-pixels within the image-domain results of the 
respective one or more discontinuities, the sub-pixels 
being transferred within their respective image-domain 
results from one side to the other of said location. 

The de-convolution process may be carried out using 
least-squares running filtering. 

An imaging method and system according to the present 
invention will now be described, by way of example, with 
reference to the accompanying drawings, in which: 

Figure 1 illustrates schematically the method and system 
of the invention; 
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Figure 2 illustrates features of processing performed in 
the method and system of Figure 1; 

Figure 3 shows results achieved from use of the method 
and system of Figure 1; 

Figure 4 shows to an enlarged scale a section of the 
contour of an image profile depicted in Figure 2; 

Figures 5 and 6 are a plan view and sectional end- 
elevation of a couch-top used in the method and system of 
Figure 1; and 

Figure 7 provides illustrates of a convolution operation, 
as a basis for a mathematical model of de-convolution 
processing in accordance with the method and system of 
Figure 1. 

The method and system to be described with reference to 
Figure 1 utilise MR scanning for medical diagnostic and 
treatment-planning purposes. In principle and in the 
general techniques described, the method and system of 
the invention can be used in other applications of MR 
scanning and also in circumstances where other scanning 
techniques are utilised. Furthermore, although both 
structure and function are represented by discrete 
x boxes' 1 to 19 in Figure 1, the method and system are to 
a substantial extent manifest in programmed digital data- 
processing operations. 

Referring to Figure 1, data derived in accordance with 
conventional operation of an MR scanner 1 is processed 
for imaging purposes within a processor 2 . The output of 
the processor 2 is used to provide a display 3, and from 
this is subject to post-imaging processing 4. The post- 
imaging processing 4 includes the facility for selecting 
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a region of the display 3 for more-detailed and closer 
inspection. 

To the extent the imaging method and system of Figure 1 
have so far been described, they are conventional, and it 
is in further processing 5 of the image data of the 
selected region of interest obtained by the post-imaging 
processing 4 that a step forward from what is already 
known is achieved* More particularly, the further 
processing 5 is operative to define more clearly the true 
edges or boundaries of image contour (s) in the selected 
region of interest, and to enhance the accuracy of the 
imaging of those contours. 

The definition and accuracy of transfer of features from 
the object domain scanned by the scanner 1, to the image 
domain manifest in the post-imaging processing 4, is 
limited by many factors. The limitations arise from 
within the scanner 1 itself (in particular the finite 
size of the energy source) , within the processing 
performed by the processor 2, and within the display 3; 
limitations arise inherently from, for example, the data 
sampling frequency and display density used, and also 
from the filter-function of the software involved. More 
particularly, and referring to Figure 2, an infinitely 
fine delta function 20 in the object domain is not 
faithfully reproduced in the image domain. Instead, the 
transfer as represented by the arrow 21 results in a 
point-spread function (PSF) or smeared-out image 22 in 
the image domain. Similarly, an infinitely sharp edge- 
response function (ERF) or step 2 3 in the object domain 
becomes through the transfer represented by arrow 24, a 
smeared-out transition 2 5 of spread represented by 
dimension arrows 26, in the image domain. When two image 
ERFs are close to one another, the smeared-out effects 
run into each other. The consequent deterioration of the 
spatial resolution is often monitored by the percentage 
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modulation transfer which is given by the ratio, 
expressed as a percentage, of the amplitude of the 
modulation in the image domain to that in the object 
domain* 

5 

The smearing effect becomes more intense as adjacent ERFs 
of discontinuities or contrast profiles get closer to 
each other (or as the spatial frequency of the modulation 
becomes higher); this also causes loss of profile height. 
10 The inherent loss of the spatial resolution (that is, the 
part that is indicated by the smeared-out effect on the 
corner edge of the ERF) cannot be restored or partially 
restored even by re-scanning the image with an ultra high 
resolution digital scanner system. 

15 

The further processing 5 is operative in accordance with 
the invention to provide accurate image-edge definition 
and location, and to improve spatial resolution in the 
imaging. More especially, in the context of Figure 2, 

20 the edge position corresponding to the discontinuity or 
step 23 of the object ERF is pin-pointed in the image 
domain from the mid-point of the full-width half-maximum 
(FWHM) of the image PSF; the pin-pointing may be to 
sub-pixel accuracy for the image ERF. In a practical 

25 implementation, low-contrast filtering, % area' filtering, 
and sub-pixel sampling may be used to remove the 
* spurious' edges and other features of the single-pixel 
modulation. The discontinuity or step 2 3 of the ERF is 
then restored within the image domain by removing the 

30 sub-pixel values from outside the optimum edge position 
to compensate for those within. It is to be noted that 
the sub-pixels then become pixels in display, and that 
the enhancement is equivalent to the performance of an 
extra high resolution image transfer system. 



As represented in Figure 2 by the arrows 2 7 and 28, the 
image ERF 25 of an infinitely-sharp ERF 23 can be 
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produced by convolution of the image PSF 2 2 with the 
object ERF 23. In accordance with the present invention, 
de-convolution of the image ERF 2 5 represented by the 
arrow 29 reproduces the image PSF 22 in a de-convolution 
5 space as image PSF 22'. The image ERF 2 5 is superimposed 
on the PSF 22' within this space as image ERF 25', and 
the optimum edge-position 3 0 is derived from the mid- 
point of the FWHM of the image PSF 22 ', and is pin- 
pointed to sub-pixel accuracy. 

10 

For one-dimensional cases, the operation in accordance 
with the invention is relatively simple, as only either 
the x- or the y-profile, that is to say a line spread 
function LSF is involved. But for two-dimensional 
15 operations, both the x- and y-profiles, and if necessary, 
the xy-diagonal profiles to eliminate any possible 
streakings in the image, may be used; in this case, a 
proper weighting scheme will be required to re-construct 
the image . 

20 

Once the original sharp-edge feature represented by the 
object ERF 23 is pin-pointed at the position 30 within 
further processing 5, that feature may be restored by 
additional re-processing 6 (Figure 1) . In re-processing 

25 6, the sub-pixel values occurring v outside' the optimum 
edge-position 30 are transferred to compensate those 
x within'. This is illustrated in Figure 2 by arrow 31 
transferring sub-pixel blocks 32 from after point 30 in 
the image ERF 25', to before it. The re-construction of 

30 image ERF 25' into image ERF 3 3 conforming closely in 
configuration to object ERF 2 3 is represented by arrow 
34. Image ERF 33 is displayed in enlarged form in 
display 7 (Figure 1) . 

35 These techniques enable substantial recovery of the loss 
of spatial resolution in the imaging, without the 
trade-off loss of other properties such as image noise. 
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Furthermore, the enhancement of spatial resolution in 
display 7 reproduces the region of interest selected from 
display 3, without blurring (or step) effects at the 
profile edge. 

5 

Figure 3 is illustrative of some of the low-contrast 
results provided in practice from MR scanning of a pig's 
brain in fluid. Curve A is the image ERF produced, 
whereas curve B is the line spread function (LSF) 

10 resulting from de-convolution of curve A carried out in 
processing 5. The optimum edge-position is established 
from the mid-point C of the FWHM of curve B, and the 
additional re-processing 6 is operative by means of sub- 
pixel transfer, to re-construct curve A to conforms 

15 substantially to the edge-feature from which it 
originated in the display 7. 

It is to be noted that whereas curve A is stepped, curve 
B is nonetheless smooth and that mid-point C is located 
20 to sub-pixel accuracy. Furthermore curve B indicates a 
sensitivity of more than 8:1 between the profile-height 
and background noise. 

Referring again to Figure 2, the complete profile 35 of 
25 an image within the selected area of interest of display 
3, is built up as indicated by arrow 36, from the edge- 
position data derived within processing 5. This data 
identifies the location of the point 30, together with 
the locations of all corresponding points derived from 
30 sampling the multiple x- or y-profiles of the selected 
area of interest. The build up and display of these 
points from the data takes place in display 8 so that a 
substantially true contour 37 for the profile 3 5 is 
defined. The sharpness of the true contour 3 7 is in 
3 5 contrast to the smeared contour that without de- 
convolution would have been obtained by virtue of the 
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spread (represented by the arrows 26) of the relevant 
image ERFs 25. 

A small portion of the contour 37 is shown enlarged in 
5 Figure 4 and is defined as a best-fit line between 

optimum edge-positions derived respectively from x- and 
y-profiles; the x-profile positions are indicated by 
black dots and the y-profile positions by white dots. 
The closeness to one another of corresponding positions 

10 in the x- and y-profiles is indicative of the accuracy to 
sub-pixel level achieved. The smear-out that would have 
been manifest in the image-profile contour if the de- 
convolution technique were not used, would have extended 
throughout the space bounded by the dashed lines 3 8 and 

15 39; these boundaries are indicated in dashed line in the 
representation of profile 35 in Figure 2. 

The accurate definition of the image contour 37 derived 
in the display 8 allows correspondingly accurate 

20 determination in calculation 9 of the area within that 
contour; the volume involved can also be derived from 
successive slices scanned. The determination of area and 
volume is especially useful for diagnostic and accurate 
assessment 10 of the size of a tumour or lesion before 

25 and after treatment. It is similarly useful for 
assessment of arterial dimensions in angiography. 

Moreover, the accurate definition of the image contour 37 
derived in the display 8, is particularly useful for 

30 segmenting anatomical structures for diagnostic and 
treatment planning 11. Furthermore, the ratio of 
intensities of two scans are derived by processing 12 to 
derive values of relaxation times T { and T 2/ and of proton 
density. The values are then represented in display 13 

35 within the boundary of the image contour, utilising 

standardisation data derived from a couch-top 14 used 
within the scanner l. The couch-top 14, which also 
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provides landmarks for determining position coordinates, 
has the form shown in Figures 5 and 6, and will now be 
described. 

5 Referring to Figures 5 and 6, the couch-top 40, which is 
of polystyrene bubble foam, has the form of a flat slab, 
and is supported on a curved foam base 41 that fits to 
the scanner table of the MR installation. Two zig-zag 
runs of tubing 4 2 are embedded within the top 4 0 to 
10 extend along the length of the couch either side of 
centrally-embedded straight tubing 43. 

The tubing 42 of each zig-zag run is of double-bore 
rectangular cross-section, whereas the tubing 43 is of 

15 single-bore cross-section. The five bores defined by the 
array of tubing 42 and 43 may be filled respectively with 
the five MR solutions S 0 to S 4 of Table I, for 
standardisation and calibration purposes. The four 
MnCl 2 .4H 2 0 solutions, S x to S 4 , cover the full range of 

20 values of T 2 and T 2 for anatomical tissues, and the fifth 
solution, S 0 , of CuS0 4 .5H 2 O, is nominally equivalent to 
"loosely bound water 11 . 



Table I 



Solution 




T, 

at 0.5T 


T 2 

at 0.5T 


S 0 


1.2 5 g/1 CuS0 4 .5H 2 O 


2 00 ms 


2 00 ms 


s, 


3.41X10 16 Mn +2 ions/ml 


840 ms 


300 ms 


s 2 


1.15X10 17 Mn +2 ions/ml 


44 0 ms 


120 ms 


s 3 


2.30X10 17 Mn +2 ions/ml 


250 ms 


60 ms 


s 4 


4.37X10 17 Mn +2 ions/ml 


150 ms 


3 0 ms 
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The tissue types revealed by the T, , T 2 and proton density 
values in display 13 are determined by processing 15 from 
look-up tables, and tissue densities are assigned within 
the image-contour boundaries in display 16. The images 
5 of display 16 are furthermore corrected empirically for 

geometry distortion in accordance with data supplied from 
memory 17, The data stored in memory 17 is derived using 
a set of drum phantoms of the spider-web form, and 
correction for geometry distortion is realistically 
10 effective principally because of the accuracy and spatial 
resolution with which image-contours are defined. 

The tissue types assigned to the corrected images are 
utilised in display 16 through the Bulk Heterogeneity 

15 Correction method described by Richard A, Geise et al , 

Radiology, 124:133-141, July, 1977, to establish for each 
image a normalised tissue density value; the up-datable 
look-up table for this is stored in memory 18. 
Accordingly, the display 16 when used at step 19 in 

20 conjunction with the positional datums derived from the 
couch 14, has all the tissue contours accurately mapped 
out with their respective tissue densities and locations. 
This establishes an accurate and readily-usable, stand- 
alone basis for diagnosis and treatment planning, and 

25 enables a true three-dimensional assessment and plan to 
be made when both orthogonal and oblique MR images are 
involved. 

Although the method and system of the invention have been 
30 described above in the medical context they are 

applicable more widely than this, for example, in 
engineering, in physical science and in the field of 
instrumentation generally. Moreover, the method and 
system is not limited to MR imaging, but may be utilised 
35 where other forms of imaging are involved. The steps and 
structure represented in Figure 1 by x boxes' 1 to 11 are 
just as applicable to computer assisted tomography (CT) , 
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as they are to MR imaging. Other forms of imaging to 
which the invention is applicable include X-ray 
radiography, film- or print-image transformation to 
digital form, digital X-ray f luorography , ultra-sound 
5 imaging, nuclear medicine, positron emission tomography 
(PET) and other camera or imaging. The technique is 
particularly suitable for use in X-ray digital 
f luorography , in which small structures under study are 
highlighted by injection of contrast liquids; the small 
10 structures may also be isolated from surrounding 
interfering effects by using an image subtraction 
technique. 

The inherent resolutions of X-ray radiography, 
15 ultra-sound imaging, nuclear medicine, and PET scanning 
are relatively low, and some are used for real-time 
study. Only the individual still frame or hard-copy 
images may be re-processed. 

In the context of engineering, physical science and the 
field of instrumentation, the invention is applicable to 
one-dimensional imaging as used, for example,, in regard 
to bar-code patterns, the spectrum of DNA analysis, iris 
patterns of eyes (for example, for identification 
purposes in commercial banking) , finger-print 
identification, and emission spectroscopy. The invention 
is also applicable to two-dimensional imaging, for 
example, in relation to images obtained by satellite or 
pattern recognition, or from a surveillance camera or 
during laboratory experimentation. As a general matter, 
the invention is applicable where there needs to be 
accurate determination of the edge position in an image 
vejrsus the true object-edge position, for the purpose, 
for example, of measurement of the positional 
displacement between object and image, distortion 
correction and manufacturing control. 



25 



30 
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As a further example of application of the present 
invention, a method and system that uses CT and MR 
imaging in conjunction with one another, will now be 
described . 

5 

The major contribution to the magnetic-resonance (MR) 
signal comes from the abundant protons content of water 
molecules and protein. It is a quantum process at the 
Larmor frequency according to the magnetic field in use. 

10 The *T,-weighted' and x T 2 -weighted' MR signals from 

protons provide contrast numbers that are relative in 
scale, whereas in CT, the X-ray absorption is a 
polychromatic attenuation process affected by the 
electron densities of all the atoms presented within the 

15 X-ray beam. There is no equation to correlate the CT 

number (or the linear attenuation coefficient, electron 
density, or tissue density) with the MR-contrast numbers; 
no direct calibration between the two types of signal is 
possible. This lack of correlation is confirmed by 

20 consideration of bone and air which are at opposite ends 
of the CT contrast (absolute) scale using water as the 
base-line reference, but which are at the same end of the 
MR-image contrast (relative) scale owing to their common 
low proton-population. 

25 

The lack of correlation between the CT and MR signals 
acts against their use in combination for imaging 
purposes, but the present invention provides a method and 
system by which the advantages of each may be utilised to 
3 0 improve image resolution and contrast information. 

In the latter regard, CT provides a high spatial 
resolution but only in regard to view normal to the 
transverse slice. Resolution for all re-constructed 
35 non-transverse planes is poor owing to the need to use 

elongate voxels to improve signal-to-noise ratio. Also, 
the partial-volume effect of using elongate voxels may 
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give rise to detection errors at the thin edge of a 
contrast profile of a lesion. MR, on the other hand, can 
give the same high degree of spatial resolution viewed in 
the normal direction to any image-slice plane, and can 
5 also provide isotropic resolution with cubic-voxel volume 
imaging. 

To this end, multiple-slice transverse CT scans are 
collected across a section of the volume of interest in a 

10 patient or other subject. Corresponding multiple-slice 
transverse MR scans of the same volume are also 
collected. The slice thickness of the latter scans may 
be one half, or smaller, of the thickness of the CT 
slices, and may be collected two-dimensionally or three- 

15 dimensionally . The patient or other object scanned is 

constrained throughout on a couch that provides landmarks 
with respect to a coordinate reference arrangement on the 
couch-top. This is to ensure the reproducibility of, for 
example, anatomical positions and features to the first 

20 order accuracy for the corresponding CT and MR scans, and 
possibly for radiation treatment to be made. 

The respective transverse planes of the CT and MR images 
are processed individually in the method and are matched 

25 with one another in a de-convoluted space for the CT and 
MR images. The two sets of de-convoluted maps are then 
merged together to a second order of accuracy in order 
that the CT numbers may be transferred over to replace 
the corresponding MR contrast numbers . Once this has 

30 been achieved, non-transverse (or oblique) planes can be 
obtained from the two-dimensional MR images, or from the 
re-arrangement of the corresponding voxels of the three- 
dimensional volume images; where two-dimensional MR is 
used, a further step of contrast transformation may be 

3 5 required. 
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The respective transverse planes of CT and MR images are 
processed individually by using x boundary' or 
x finger-print ' matching techniques in a de-convolution 
space for the CT & MR images. In these transverse image 
5 planes, in particular when used medically, the 

skin-contour features along the sides of the patient may 
be best used for second-order alignment and matching 
purposes, as they are less affected by patient-movements. 
Transitional error may be readily corrected with respect 

10 to the coordinate positions of a rectangular tubing 
system embedded in the couch-top (for example that 
described above with reference to Figures 5 and 6) . The 
processed data may then be used for a x diagnostic and 
statistics software package' of CT image versus MR image 

15 for their exactly corresponding transverse slice (s), and 
an associated statistical package' for accurate 
computation of the x true' area, and then the 'true' 
volume, of a lesion or tissue profile or contour. 

20 The two sets of de-convoluted maps may also be merged 

together to a second order of accuracy in order that the 
CT numbers may be transferred over to replace the 
corresponding MR contrast numbers. Once this has been 
done, non-transverse (or oblique) planes are obtained 

25 from the two-dimensional image or from the re-arrangement 
of the corresponding voxels of the three-dimensional 
volume; in the two-dimensional MR case, a further step of 
contrast transformation is required. The transferred 
contrast data may then be used in a three-dimensional 

30 radiotherapy treatment planning software package for an 

in-plane, oblique-image pseudo-three-dimensional approach 
using MR images. 

Software required for de-convolution processing of image 
3 5 data according to the invention may be implemented in 

conjunction with a least-squares curve-fitting method. A 
mathematical model of the method, from which the required 
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software can be readily developed, will now be given in 
relation to Figure 7 which indicates the convolution of 
an LSF L(x) with a step function SF having a step-down 
edge SDE, to produce an image ERF E(x) . Considering the 
values of E(x) resulting from the convolution at points 
x=l to x=n of SF: 



15 



at x=l : 



£l,U) Ax = E(x) x 



10 at x=2 



^L(x) Ax - ]Tl,(x) Ax = E ^2 



at x=3 



L <x) Ax - L {x) Ax = E{x) 3 



a-2 



until at x=n: a 

£l(x) Ax - £l,(x) Ax = £(x) fl 



For de-convolution, that is to say, for the ERF E(x) to 
manifest the LSF L(x) : 



E(X) 1 - E(X) 2 = J2 L{X) AX = LU) (a-D-a Ax 
a-l 

{E(X) X - £(x) 2 J/Ax = L(X) (a.D-a 

thus : 

{£ , (X) 2 - £<X) 3 )/AX = L(X) (a _ 2 )-<a-l) 

leading to 

te(x> - i?(x) n )/Ax = L(x) 



20 
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Thus, the same shape of LSF is recovered by the 
de-convolution process independently of the sense of the 
ERF; it does not matter whether a roll-up or a roll-down 
ERF of a positive or a negative contrast profile, is 
5 involved. This is an important property as a positive 
contrast contour will have roll-up (from low to high) 
ERFs at both ends whereas a negative contrast contour 
will have roll-down (from high to low) ERFs at both ends. 

10 In the practical implementation, the LSF is derived from 
the ERF by de-convolution using a running filter. This 
enhances the accuracy of the method in overcoming the 
problem of noise that affects the digitised pixel values 
of the image. Use is made of a least-squares fitting of 

15 a set of data points that advances along the whole length 
of the function from one point to another. 



20 



Assuming that y = a 0 + a,x + a 2 x 2 represents the ERF curve 
a five-point or seven-point fit is used, and the normal 
equation becomes: 



Ei Ex.. Exj 
Ex_£ Ex? Exj 



aO 
al 
a2 



Ex iyi 



25 



where: i is 1+n, ... 5+n; ... until m-5, ... m 
or 

i is 1+n, . . . 7+n; . . . until in— 7 , . . . m 
n is 0, 1 , 2 , ... ; and 

1, ... m is the span of the ERF profile, 



The solution may be derived from either: 



aO 




Ei 


Ex 


Ex 2 


-1 


Ey 


al 




Ex 


Ex 2 


Ex 3 




Exy 


a2 




Ex 2 


Ex 3 


Ex 4 




Ex 2 y 
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or: 



aO 



Ey Ex Ex 2 
Exy Ex 2 Ex 3 
Ex 2 y Ex 3 Ex 4 



El Ey Ex 
Ex Exy Ex 3 
Ex 2 Ex 2 y Ex 4 



a2 



El Ex Ey 
Ex Ex 2 Exy 
Ex 2 Ex 3 Ex 2 y 



El Ex Ex 2 
Ex Ex 2 Ex 3 



-X" 



.X 4 



For both of these equations to be valid :- 



El Ex Ex 2 
Ex Ex 2 Ex 3 
Ex 2 Ex 3 Ex 4 



* 0 



5 The gradient at dy/dx (3+n) or dy/dx (4+n) can then be derived 
and plotted against x for the LSF profile. 

The graph of dy/dx against x gives the LSF profile. The 
peak of this profile is located centrally of the mid- 
10 points of the ascending and descending limbs of the 

graph. These points define the extremes of the full- 
width-half -maximum (FWHM) of the profile and the mid- 
point of this is determined with an accuracy of sub-pixel 
level owing to the * average' effect. 

15 

The point spread function (PSF) is the two dimensional 
profile which may be derived, in practice, from the two 
corresponding LSFs orthogonal to one another within an 
image plane. The peak position of the PSF profile is, 
20 therefore, from the x mean' or x cross-over' of the two 
peaks or the two LSF profiles. The PSF is obtained in 
practice from two orthogonal axes in a two-dimensional 
plane. 



25 The generation of the LSF (or PSF) is, after phase- 
reversal correction, independent of the roll-up or 
roll-down nature of ERFs at the edges of the contrast 
contour. In other words, it is independent of the sense 
and the absolute value of the contrast numbers within the 
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contour. The peak position of the LSF (or PSF) is the 
central half-way (50%) point of the roll-up or roll-down 
ERF which is the optimum position for true-edge 
definition. 
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Claims : 

1. A method of imaging wherein a de-convolution process 
is applied to the image-domain results of an object-scan 
to derive therefrom the respective point- or line-spread 
function effective in the object- to image-domain 
transfer of one or more object-discontinuities, and to 
derive from said function the location in the image 
domain of the respective discontinuity. 

2* A method according to Claim 1 wherein the location 
of the respective discontinuity is derived from the mid- 
point of the full-width half -maximum of said function. 

3. A method according to Claim 1 or Claim 2 wherein 
said function is correlated with the image-domain results 
of said transfer for enhancement of spatial resolution of 
the imaging of the one or more discontinuities. 

4. A method according to Claim 3 wherein the 
enhancement of spatial resolution of the imaging of the 
one or more discontinuities involves transfer of sub- 
pixels within the image-domain results of the respective 
one or more discontinuities, the sub-pixels being 
transferred within their respective image-domain results 
from one side to the other of said location. 

5. A method according to any one of Claims 1 to 4 
wherein the de-convolution process is carried out using 
least-squares running filtering. 

6. A method according to any one of Claims 1 to 5 
wherein an edge-contour of the object is defined in the 
image domain using de-convolution processing as 
aforesaid. 
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7. A method according to Claim 6 wherein the area 
and/or volume of the object-image within the edge-contour 
is determined. 

8- A method according to Claim 6 or Claim 7 wherein the 
object-scan is a magnetic resonance (MR) scan, values of 
relaxation times T, and T 2 are derived for the object- 
image within said contour, and these values are used to 
identify from stored data, types of tissue or other 
material involved in the scanned object. 

9. A method according to Claim 8 wherein density values 
for the identified tissue or other material types are 
derived from further stored data. 

10. A method of imaging according to any one of Claims 1 
to 9 wherein corresponding computed tomography (CT) and 
magnetic resonance (MR) scans of the same part of an 
object are derived, the scans are related to one another 
for correlation of one to the other positionally with 
respect to said part using the de-convolution process, 
and imaging of said part of the object is provided in 
accordance with the MR scan as modified spatially in 
dependence upon the CT contrast numbers applicable to the 
corresponding, correlated positions of the CT scan. 

11. A method according to any one of Claims 8 to 10 
wherein geometric correction is applied to the imaging 
derived from the MR scan, in accordance with stored data. 

12. An imaging system comprising means for performing a 
de-convolution process on the image-domain results of an 
object-scan to derive therefrom the respective point- or 
line-spread function effective in the object- to image- 
domain transfer of one or more object-discontinuities, 
and means to derive from said function the location in 
the image domain of the respective discontinuity. 
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13. A system according to Claim 12 wherein the location 
of the respective discontinuity is derived from the mid- 
point of the full-width half-maximum of said function. 

14. A system according to Claim 12 or Claim 13 wherein 
said function is correlated with the image-domain results 
of said transfer for enhancement of spatial resolution of 
the imaging of the one or more discontinuities. 

15. A method according to Claim 14 wherein the 
enhancement of spatial resolution of the imaging of the 
one or more discontinuities involves transfer of sub- 
pixels within the image-domain results of the respective 
one or more discontinuities, the sub-pixels being 
transferred within their respective image-domain results 
from one side to the other of said location. 

16. A system according to any one of Claims 12 to 15 
wherein the de-convolution process is carried out using 
least-squares running filtering. 
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I. Basis of the report 

1 . This report has been drawn on the basis of (substitute sheets which have been furnished to the receiving Office in 
response to an invitation under Article 14 are referred to in this report as "originally filed" and are not annexed to 
the report since they do not contain amendments (Rules 70. 16 and 70.17).): 

Description, pages: 

3,4,8-1 8 as originally filed 

1 ,2,5-7 as received on 22/1 1/2000 with letter of 1 8/1 1/2000 

Claims, No.: 

1-16 as received on 22/1 1 /2000 with letter of 1 8/1 1 /2000 

Drawings, sheets: 

1 ,2,4-6 as originally filed 

3 as received on 22/1 1/2000 with letter of 18/1 1/2000 

2. With regard to the language, all the elements marked above were available or furnished to this Authority in the 
language in which the international application was filed, unless otherwise indicated under this item. 

These elements were available or furnished to this Authority in the following language: , which is: 

□ the language of a translation furnished for the purposes of the international search (under Rule 23.1 (b)). 

□ the language of publication of the international application (under Rule 48.3(b)). 

□ the language of a translation furnished for the purposes of international preliminary examination (under Rule 
55.2 and/or 55.3). 

3. With regard to any nucleotide and/or amino acid sequence disclosed in the international application, the 
international preliminary examination was carried out on the basis of the sequence listing: 

□ contained in the international application in written form. 

□ filed together with the international application in computer readable form. 

□ furnished subsequently to this Authority in written form. 

□ furnished subsequently to this Authority in computer readable form. 

□ The statement that the subsequently furnished written sequence listing does not go beyond the disclosure in 
the international application as filed has been furnished. 

□ The statement that the information recorded in computer readable form is identical to the written sequence 
listing has been furnished. 
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4. The amendments have resulted in the cancellation of: 

□ the description, pages: 

□ the claims, Nos.: 

□ the drawings, sheets: 

5. □ This report has been established as if (some of) the amendments had not been made, since they have been 

considered to go beyond the disclosure as filed (Rule 70.2(c)): 

(Any replacement sheet containing such amendments must be referred to under item 1 and annexed to this 
report.) 



6. Additional observations, if necessary: 



V. Reasoned statement under Article 35(2) with regard to novelty, inventive step or industrial applicability; 
citations and explanations supporting such statement 

1. Statement 



Novelty (N) 


Yes: 


Claims 


4,15 




No: 


Claims 


1-3, 6-7, 12-14, 16 


Inventive step (IS) 


Yes: 


Claims 


4,15 




No: 


Claims 


5, 8-1 1 


Industrial applicability (IA) 


Yes: 


Claims 


1-16 




No: 


Claims 





2. Citations and explanations 
see separate sheet 



VIII. Certain observations on the international application 

The following observations on the clarity of the claims, description, and drawings or on the question whether the 
claims are fully supported by the description, are made: 
see separate sheet 
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Re. Sections V, VIII (Novelty, Inventive step, Observations) 

1 . The following documents are referred to in this report:- 

D1 - D3: The documents of the International Search Report in the order as 
they are listed therein. 



2. Claim 1 is not clear in various respects (Art. 6 PCT):- 

2.1 Claim 1 specifies in the first lines that "a de-convolution process is applied to 
image-domain results of an object-scan to derive therefrom the respective point- 
or line-spread function of ... object discontinuities". This phrase has not been 
understood. While the first part of the phrase can be well understood as to 
address, in generic terms, a de-convolution process applied to image-domain 
results of an object-scan, it is not at all clear, how from a scan of a real-world 
object a point- or line-spread function can be obtained, as suggested in the 
second part of the criticized phrase. 

Conventionally, the point- or line-spread function of an imaging system is 
determined prior to the actual object-scans of interest. The point- or line-spread 
function is then used, in a de-convolution step, to correct the object scans (see 
e.g. D1, page 154, or D2, abstract and page 153). Since, from pages 4-6 and from 
Fig.2 of the present application it appears that this basic de-convolution technique 
is also used in the context of the claimed invention, current claim 1 seems 
misleading in this respect. 

2.2 It is not clear what exactly is meant with the last feature "sub-pixel sampling". The 
feature is therefore interpreted to merely mean that the whole de-convolution 
process is carried out at sub-pixel accuracy. 

2.3 In the light of the description and on basis of general knowledge (see in particular 
D1 and D2) claim 1 is therefore interpreted as to recite: 

A method of imaging in which a de-convolution process is applied to image- 
domain results of an object-scan, using an imaging system, the method 
comprising the steps of> 
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j (a) providing a point- or line-spread function (22) of the imaging system; 
- £ (b) ^canning an object having a discontinuity, thereby obtaining an image of the 
object; 

(c) deriving from the point- or line-spread function (22) the location of the 
discontinuity in the image, 

(d) wherein the de-convolution process is carried out at sub-pixel accuracy. 



3. Novelty and Inventive step. 

3.1 The feature (c) of claim 1 ("to derive from said function the location ...") is broad 
and clearly reads onto the conventional deconvolution method in which the scan 
image is de-convolved with the LSF or PSF of the imaging system; this usually 
resulting in an edge-enhanced image from which "object discontinuities", i.e. 
edges, can be "derived", exactly as it is worded in the claim, see e.g. D1: Fig. 3, or 
D2: abstract and page 1 53. 

The last feature "sub-pixel accuracy" is known from D1, page 154, right hand 
column, lines 8-1 1 . 

It thus appears that claim 1 lacks novelty over D1 (Art.33(2) PCT). 

3.2 Document D3 discloses edge feature extraction at sub-pixel accuracy (see 
abstract; page 344, col.1), and a combination of D2 and D3 would seem to render 
claim 1 obvious. 

It is moreover noted that, further to D3, various techniques of "sub-pixel sampling" 
of some kind are known in the art to improve the perceived image sharpness, and 
it is not apparent in what respect the broad claimed feature "sub-pixel sampling" 
would go beyond these conventional techniques. 



4. Dependent claims. 

- Claim 2: The feature "mid-point of the full-width half-maximum" of the LSF 
seems to be known from D1: page 155, col. 2 and D2: page 157, col.1 and table II. 
Should there be any dispute on the disclosure of D1 and D2 in this respect, then 
the feature at issue would at least lack an inventive step. 

- Claim 3: The limitation of the claim with respect to claim 1 and conventional 
deconvolution methods is not apparent. 
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- Claim 4: The notion of transfer of sub-pixels is not known from the prior art on 
file. The solution provided by claim 4 seems to enable substantial recovery of the 
loss of spatial resolution in the imaging, without the trade-off loss of other 
properties such as image noise which is often a problem in edge enhancement. 
While D3 is also concerned with edge feature extraction at sub-pixel accuracy, D3 
uses another approach and classifies an edge with respect to edge profiles (D3: 
abstract; page 344, col.1). 

- Claim 5: Least-square filters are used in various contexts (see D2: abstract; D3: 
abstract) and the claim is too vague to provide any new effect. 

Claims 6, 7: see D1, Fig. 3. 

Claims 8-1 1 : The technical problems solved by subject-matter of these claims 
seem unrelated with that of the preceding claims in that it concerns quite different 
aspects of an MR or CT system. The features of the claims as such seem to be 
conventional in the art, and the claims therefore seem to lack an inventive step. 



5. The above objections apply to the system claims accordingly. 

6. Briefly, it seems that the contribution of the application over the prior art seems to 
reside in the particular way in which unsharp edges are restored, i.e. in the 
algorithm that uses transfer of sub-pixels, as specified in dependent claim 4. 

***** 
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modulation transfer which is given by the ratio, 
expressed as a percentage, of the amplitude of the 
modulation in the image domain to that in the object 
domain. 

The smearing effect becomes more intense as adjacent ERFs 
of discontinuities or contrast profiles get closer to 
each other (or as the spatial frequency of the modulation 
becomes higher); this also causes loss of profile height. 
The inherent loss of the spatial resolution (that is, the 
part that is indicated by the smeared-out effect on the 
corner edge of the ERF) cannot be restored or partially 
restored even by re-scanning the image with an ultra high 
resolution digital scanner system. 

The further processing 5 is operative in accordance with 
the invention to provide accurate image-edge definition 
and location, and to improve spatial resolution in the 
imaging. More especially, in the context of Figure 2, 
the edge position corresponding to the discontinuity or 
step 2 3 of the object ERF is pin-pointed in the image 
domain from the mid-point of the full-width half-maximum 
(FWHM) of the image PSF; the pin-pointing is to sub-pixel 
accuracy for the image ERF. Low-contrast and >area' 
filtering are used to remove x spurious' edges, and sub- 
pixel sampling to detect detail to the resolution of the 
single-pixel modulation. The discontinuity or step 2 3 of 
the ERF is then restored within the image domain by 
removing the sub-pixel values from outside the optimum 
edge position to compensate for those within. It is to 
be noted that the sub-pixels then become pixels in 
display, and that the enhancement is equivalent to the 
performance of an extra high resolution image transfer 
system . 

As represented in Figure 2 by the arrows 27 and 28, the 
image ERF 25 of an infinitely-sharp, step 2 3 can be 
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produced by convolution of the image PSF 2 2 with the 
object ERF 23. In accordance with the present invention, 
de-convolution of the image ERF 2 5 using sub-pixel 
sampling, represented by the arrow 29, reproduces the 
5 image PSF 22 in a de-convolution space as image PSF 22 
The image ERF 25 is superimposed on the PSF 22' within 
this space as image ERF 25', and the optimum edge- 
position 3 0 is derived from the mid-point of the FWHM of 
the image PSF 22', and is pin-pointed to sub-pixel 
10 accuracy. 

For one-dimensional cases, the operation in accordance 
with the invention is relatively simple, as only either 
the x- or the y-profile, that is to say a line spread 
15 function LSF is involved. But for two-dimensional 

operations, both the x- and y-profiles, and if necessary, 
the xy-diagonal profiles to eliminate any possible 
streakings in the image, may be used; in this case, a 
proper weighting scheme will be required to re-construct 

2 0 the image. 

Once the original sharp-edge feature represented by the 
object ERF 23 is pin-pointed at the position 3 0 within 
further processing 5, that feature may be restored by 
25 additional re-processing 6 (Figure l) . In re-processing 
6, the sub-pixel values occurring x outside' the optimum 
edge-position 3 0 are transferred to compensate those 
^within'. This is illustrated in Figure 2 by arrow 31 
transferring sub-pixel blocks 32 from after point 30 in 

3 0 the image ERF 25', to before it. The re-construction of 

image ERF 25' into image ERF 3 3 conforming closely in 
configuration to object ERF 2 3 is represented by arrow 
34. Image ERF 3 3 is displayed in enlarged form in 
display 7 (Figure 1) . 



These techniques enable substantial recovery of the loss 
of spatial resolution in the imaging, without the 
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trade-off loss of other properties such as image noise. 
Furthermore, the enhancement of spatial resolution in 
display 7 reproduces the region of interest selected from 
display 3, without blurring- (or step) effects at the 
5 profile edge. 



Figure 3 is illustrative of some of the low-contrast 
results provided in practice from MR scanning of a pig's 
brain in fluid. Curve A is the image ERF produced, 

10 whereas curve B is the line spread function (LSF) 

resulting from de-convolution of curve A carried out in 
processing 5. The optimum edge-position is established 
from the mid-point C of the FWHM of curve B, and the 
additional re-processing 6 is operative by means of sub- 

15 pixel transfer, to re-construct curve A to conform 
substantially to the edge-feature from which it 
originated in the display 7. 

It is to be noted that whereas curve A is stepped, curve 
20 B is nonetheless smooth and that mid-point C is located 
to sub-pixel accuracy. Furthermore curve B indicates a 
sensitivity of more than 8:1 between the profile-height 
and background noise. 

25 Referring again to Figure 2, the complete profile 35 of 
an image within the selected area of interest of display 
3, is built up as indicated by arrow 36, from the edge- 
position data derived within processing 5. This data 
identifies the location of the point 30, together with 

3 0 the locations of all corresponding points derived from 
sampling the multiple x- or y-profiles of the selected 
area of interest. The build up and display of these 
points from the data takes place in display 8 so that a 
substantially true contour 37 for the- profile 35 is 

3 5 defined. The sharpness of the true contour 37 is in 
contrast to the smeared contour that without de- 
convolution would have been obtained by virtue of the 
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Claims: 

1. A method of imaging in which a de-convolution 
process is applied to the image-domain results of an 
object-scan to derive therefrom the respective point- or 
line-spread function of one or more object- 
discontinuities, and to derive from said function the 
location in the image domain of the respective 
discontinuity, wherein the de-convolution process is 
carried out using sub-pixel sampling. 

2. A method according to Claim 1 wherein the location 
of the respective discontinuity is derived from the mid- 
point of the full-width half -maximum of said function. 

3. A method according to Claim 1 or Claim 2 wherein 
said function is correlated with the image-domain results 
of said transfer for enhancement of spatial resolution of 
the imaging of the one or more discontinuities. 

4. A method according to Claim 3 wherein the 
enhancement of spatial resolution of the imaging of the 
one or more discontinuities involves transfer of sub- 
pixels within the image-domain results of the respective 
one- or more discontinuities, the sub-pixels being 
transferred within their respective image-domain results 
from one side to the other of said location for edge- 
image definition. 

5. A method according to any one of Claims 1 to 4 
wherein the de-convolution process is carried out using 
least-squares running filtering. 

6. A method according to any one of - Claims 1 to 5 
wherein an edge-contour of the object is defined in the 
image domain using de-convolution processing as 
aforesaid. 
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7. A method according to Claim 6 wherein the area 
and/or volume of the object-image within the edge-contour 
is determined. 

8. A method according to Claim 6 or Claim 7 wherein the 
object-scan is a magnetic resonance (MR) scan, values of 
relaxation times T x and T 2 are derived for the object- 
image within said contour, and these values are used to 
identify from stored data, types of tissue or other 
material involved in the scanned object, 

9. A method according to Claim 8 wherein density values 
for the identified tissue or other material types are 
derived from further stored data. 

10. A method of imaging according to any one of Claims 1 
to 9 wherein corresponding computed tomography (CT) and 
magnetic resonance (MR) scans of the same part of an 
object are derived, the scans are related to one another 
for correlation of one to the other positionally with 
respect to said part using the de-convolution process, 
and imaging of said part of the object is provided in 
accordance with the MR scan as modified spatially in 
dependence upon the CT contrast numbers applicable to the 
corresponding, correlated positions of the CT scan. 

11. A method according to any one of Claims 8 to 10 
wherein geometric correction is applied to the imaging 
derived from the MR scan, in accordance with stored data. 

12. An imaging system comprising means for performing a 
de-convolution process on the image-domain results of an 
object-scan to derive therefrom the respective point- or 
line-spread function of one or more object- 
discontinuities, and means to derive from said function 
the location in the image domain of the respective 



AMENDED SHEET 



GB 009903417 




21 

discontinuity, wherein the de-convolution process is 
carried out using sub-pixel sampling. 

13. A system according to Claim 12 wherein the location 
of the respective discontinuity is derived from the mid- 
point of the full-width half -maximum of said function. 

14. A system according to Claim 12 or Claim 13 wherein 
said function is correlated with the image-domain results 
of said transfer for enhancement of spatial resolution of 
the imaging of the one or more discontinuities. 

15. A system according to Claim 14 wherein the 
enhancement of spatial resolution of the imaging of the 
one or more discontinuities involves transfer of sub- 
pixels within the image-domain results of the respective 
one or more discontinuities, the sub-pixels being 
transferred within their respective image-domain results 
from one side to the other of said location for edge- 
image definition. 

16. A system according to any one of Claims 12 to 15 
wherein the de-convolution process is carried out using 
least-squares running filtering. 
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Imaging 



5 This invention relates to imaging and in particular to 
methods and systems for image enhancement. 

Imaging involves transfer from the object domain into the 
image domain, but owing to limiting factors such as the 

10 finite size of energy source, detector size, sampling 

frequency, display density, software filter function, and 
possibly partial-volume effects experienced with some 
imagers, an infinitely fine delta function in the object 
domain cannot be faithfully reproduced in the image 

15 domain. Instead, a smeared-out image, or point-spread 
function (PSF) , is observed. Similarly, an infinitely 
sharp edge-response function (ERF) in the object domain 
becomes a smeared-out ERF in the image domain. The 
smearing effect becomes more intense as the adjacent ERFs 

20 of discontinuities or contrast profiles get closer to 
each other . 

It is an object of the present invention to provide a 
method and system by which the above problem can be at 
25 least .partly overcome. 

According to one aspect of the present invention there is 
provided a method wherein a de-convolution process is 
applied to the image-domain results of an object-scan to 
30 derive therefrom the respective point- or line-spread 
function effective in the object- to image-domain 
transfer of one or more object-discontinuities, and to 
derive from said function the location in the image 
domain of the respective discontinuity. _ 

35 



According to another aspect of the invention there is 
provided an imaging system comprising means for 
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performing a de-convolution process on the image-domain 
results of an object-scan to derive therefrom the 
respective point- or line-spread function effective in 
the object- to image-domain transfer of one or more 
object-discontinuities, and means to derive from said 
function the location in the image domain of the 
respective discontinuity . 

The method and system of the invention enable the 
location of the respective discontinuity in the image 
domain, to be established with a high degree of accuracy. 
This is critical to image definition free of any 
substantial smearing, and to this end the location of the 
respective discontinuity may be derived to sub-pixel 
accuracy simply from the mid-point of the full-width 
half-maximum of said function. 

The said function may be correlated with the image-domain 
results of said transfer for enhancement of spatial 
resolution of the imaging of the one or more 
discontinuities. This enhancement may involve transfer 
of sub-pixels within the image-domain results of the 
respective one or more discontinuities, the sub-pixels 
being transferred within their respective image-domain 
results from one side to the other of said location. 

The de-convolution process may be carried out using 
least-squares running filtering. 

An imaging method and system according to the present 
invention will now be described, by way of example, with 
reference to the accompanying drawings, in which: 

Figure 1 illustrates schematically the method and system 
of the invention; 
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modulation transfer which is given by the ratio, 
expressed as a percentage, of the amplitude of the 
modulation in the image domain to that in the object 
domain, 

5 

The smearing effect becomes more intense as adjacent ERFs 
of discontinuities or contrast profiles get closer to 
each other (or as the spatial frequency of the modulation 
becomes higher); this also causes loss of profile height. 
10 The inherent loss of the spatial resolution (that is, the 
part that is indicated by the smeared-out effect on the 
corner edge of the ERF) cannot be restored or partially 
restored even by re-scanning the image with an ultra high 
resolution digital scanner system. 

15 

The further processing 5 is operative in accordance with 
the invention to provide accurate image-edge definition 
and location, and to improve spatial resolution in the 
imaging. More especially, in the context of Figure 2, 

20 the edge position corresponding to the discontinuity or 
step 23 of the object ERF is pin-pointed in the image 
domain from the mid-point of the full-width half-maximum 
(FWHM) of the image PSF; the pin-pointing may be to 
sub-pixel accuracy for the image ERF. In a practical 

25 implementation, low-contrast filtering, x area' filtering, 
and sub-pixel sampling may be used to remove the 
x spurious' edges and other features of the single-pixel 
modulation. The discontinuity or step 2 3 of the ERF is 
then restored within the image domain by removing the 

30 sub-pixel values from outside the optimum edge position 
to compensate for those within. It is to be noted that 
the sub-pixels then become pixels in display, and that 
the enhancement is equivalent to the performance of an 
extra high resolution image transfer system. 



As represented in Figure 2 by the arrows 2 7 and 28, the 
image ERF 25 of an infinitely-sharp ERF 2 3 can be 
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produced by convolution of the image PSF 22 with the 
object ERF 23. In accordance with the present invention, 
de-convolution of the image ERF 25 represented by the 
arrow 29 reproduces the image PSF 22 in a de-convolution 
5 space as image PSF 22'. The image ERF 25 is superimposed 
on the PSF 22' within this space as image ERF 25' , and 
the optimum edge-position 3 0 is derived from the mid- 
point of the FWHM of the image PSF 22' , and is pin- 
pointed to sub-pixel accuracy, 

10 

For one-dimensional cases, the operation in accordance 
with the invention is relatively simple, as only either 
the x- or the y-profile, that is to say a line spread 
function LSF is involved. But for two-dimensional 
15 operations, both the x- and y-profiles, and if necessary, 
the xy-diagonal profiles to eliminate any possible 
streakings in the image, may be used; in this case, a 
proper weighting scheme will be required to re-construct 
the image. 

20 

Once the original sharp-edge feature represented by the 
object ERF 23 is pin-pointed at the position 30 within 
further processing 5 , that feature may be restored by 
additional re-processing 6 (Figure 1) . In re-processing 

25 6, the sub-pixel values occurring x outside' the optimum 
edge-position 30 are transferred to compensate those 
x within / . This is illustrated in Figure 2 by arrow 31 
transferring sub-pixel blocks 32 from after point 30 in 
the image ERF 25', to before it. The re-construction of 

30 image ERF 25' into image ERF 3 3 conforming closely in 
configuration to object ERF 2 3 is represented by arrow 
34. Image ERF 33 is displayed in enlarged form in 
display 7 (Figure 1) . 

35 These techniques enable substantial recovery of the loss 
of spatial resolution in the imaging, without the 
trade-off loss of other properties such as image noise. 
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Furthermore, the enhancement of spatial resolution in 
display 7 reproduces the region of interest selected from 
display 3, without blurring (or step) effects at the 
profile edge, 

5 

Figure 3 is illustrative of some of the low-contrast 
results provided in practice from MR scanning of a pig's 
brain in fluid. Curve A is the image ERF produced, 
whereas curve B is the line spread function (LSF) 

10 resulting from de-convolution of curve A carried out in 
processing 5. The optimum edge-position is established 
from the mid-point C of the FWHM of curve B, and the 
additional re-processing 6 is operative by means of sub- 
pixel transfer, to re-construct curve A to conforms 

15 substantially to the edge-feature from which it 
originated in the display 7. 

It is to be noted that whereas curve A is stepped, curve 
B is nonetheless smooth and that mid-point C is located 

2 0 to sub-pixel accuracy. Furthermore curve B indicates a 

sensitivity of more than 8:1 between the profile-height 
and background noise. 

Referring again to Figure 2, the complete profile 35 of 
25 an image within the selected area of interest of display 
3, is built up as indicated by arrow 36, from the edge- 
position data derived within processing 5. This data 
identifies the location of the point 30, together with 
the locations of all corresponding points derived from 
30 sampling the multiple x- or y-profiles of the selected 
area of interest. The build up and display of these 
points from the data takes place in display 8 so that a 
substantially true contour 37 for the profile 35 is 
defined. The sharpness of the true corvtour 3 7 is in 

3 5 contrast to the smeared contour that without de- 

convolution would have been obtained by virtue of the 



WO 00/22573 




PCT/GB99/03417 



1.9 



Claims : 

1. A method of imaging wherein a de-convolution process 
is applied to the image-domain results of an object-scan 
to derive, therefrom the respective point- or line-spread 
function effective in the object- to image-domain 
transfer of one or more object-discontinuities, and to 
derive from said function the location in the image 
domain of the respective discontinuity. 

2. A method according to Claim 1 wherein the location 
of the respective discontinuity is derived from the mid- 
point of the full-width half-maximum of said function. 

3. A method according to Claim 1 or Claim 2 wherein 
said function is correlated with the image-domain results 
of said transfer for enhancement of spatial resolution of 
the imaging of the one or more discontinuities. 

4. A method according to Claim 3 wherein the 
enhancement of spatial resolution of the imaging of the 
one or more discontinuities involves transfer of sub- 
pixels within the image-domain results of the respective 
one or more discontinuities, the sub-pixels being 
transferred within their respective image-domain results 
from one side to the other of said location. 

5. A method according to any one of Claims 1 to 4 
wherein the de-convolution process is carried out using 
least-squares running filtering. 

6. A method according to any one of Claims 1 to 5 
wherein an edge-contour of the object is defined in the 
image domain using de-convolution processing as 
aforesaid. 
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7. A method according to Claim 6 wherein the area 
and/or volume of the object-image within the edge-contour 
is determined. 

8. A method according to Claim 6 or Claim 7 wherein the 
object-scan is a magnetic resonance (MR) scan, values of 
relaxation times T x and T 2 are derived for the object- 
image within said contour, and these values are used to 
identify from stored data, types of tissue or other 
material involved in the scanned object. 

9. A method according to Claim 8 wherein density values 
for the identified tissue or other material types are 
derived from further stored data. 

10. A method of imaging according to any one of Claims 1 
to 9 wherein corresponding computed tomography (CT) and 
magnetic resonance (MR) scans of the same part of an 
object are derived, the scans are related to one another 
for correlation of one to the other positionally with 
respect to said part using the de-convolution process, 
and imaging of said part of the object is provided in 
accordance with the MR scan as modified spatially in 
dependence upon the CT contrast numbers applicable to the 
corresponding, correlated positions of the CT scan. 

11. A method according to any one of Claims 8 to 10 
wherein geometric correction is applied to the imaging 
derived from the MR scan, in accordance with stored data. 

12. An imaging system comprising means for performing a 
de-convolution process on the image-domain results of an 
object-scan to derive therefrom the respective point- or 
line-spread function effective in the object- to image- 
domain transfer of one or more object-discontinuities, 
and means to derive from said function the location in 
the image domain of the respective discontinuity. 
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13. A system according to Claim 12 wherein the location 
of the respective discontinuity is derived from the mid- 
point of the full-width half-maximum of said function. 

14. A system according to Claim 12 or Claim 13 wherein 
said function is correlated with the image-domain results 
of said transfer for enhancement of spatial resolution of 
the imaging of the one or more discontinuities. 

15. A method according to Claim 14 wherein the 
enhancement of spatial resolution of the imaging of the 
one or more discontinuities involves transfer of sub- 
pixels within the image-domain results of the respective 
one or more discontinuities, the sub-pixels being 
transferred within their respective image-domain results 
from one side to the other of said location. 

16. A system according to any one of Claims 12 to 15 
wherein the de-convolution process is carried out using 
least-squares running filtering. 
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